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Motion of cantilever near sample surfaces exhibits additional friction even before two bodies 
come into mechanical contact. Called non-contact friction (NCF), this friction is of great practical 
importance to the ultrasensitive force detection measurements. Observed large NCF of a micron- 
scale cantilever found anomalously large damping that exceeds theoretical predictions by 8-11 orders 
of magnitude. This finding points to contribution beyond fluctuating electromagnetic fields within 
van der Waals approach. Recent experiments reported by Saitoh et al. (Phys. Rev. Lett. 105, 
236103 (2010)) also found nontrivial distance dependence of NCF. Motivated by these observations, 
we propose a mechanism based on the coupling of cantilever to the relaxation dynamics of surface 
defects. We assume that the surface defects couple to the cantilever tip via spin-spin coupling and 
their spin relaxation dynamics gives rise to the backaction terms and modifies both the friction 
coefficient and the spring constant. We explain the magnitude, as well as the distance dependence 
of the friction due to these backaction terms. Reasonable agreement is found with the experiments. 



Introduction: Friction is one of the most widely per- 
ceived but least understood phenomena in nature. Fric- 
tion is ubiquitously seen around us at the macroscale. 
One might expect however that new insights will be 
gained by investigating friction mechanisms at small scale 
and short distances. Indeed, recent advances in nan- 
otechnology have enabled the study of friction on the 
nanoscaleP where a novel form of friction has been dis- 
covered, namely the non- contact friction (NCF) I P3 This 
kind of friction occurs when two objects are in close prox- 
imity but not in physical contact. NCF is of great practi- 
cal importance for the modern development of ultrasen- 
sitive force detection devices The precision of these 
measurements may be ultimately limited by the effects 
of NCF. 

The origin of NCF is still under debate. One pro- 
posal was that it is the friction resulting from Ohmic 
losses mediated by fluctuating electromagnetic fields. It 
turns out be 8-11 orders of ma gnitude smaller than that 
observed in experiments! 2 * 3 * 7 *^ Several alternative mech- 
anisms have been proposed in recent years (see^ and 
references therein), but the problem remains unsolved. 




FIG. 1: (Color online) Illustration of proposed mechanism 
that gives rise to the non-contact friction between the can- 
tilever tip and the sample surface: randomly distributed de- 
fect spins on the sample surface interacting with the spins 
residing on the tip. 



More recently, systematic studies of NCF between a can- 
tilever tip and the sample surface have been performed 
on metallic, insulating and superconducting materials at 
different temperatures using hard cantilevers^ (as op- 
posed to soft cantelevers results^ ) . These studies found 
that at low temperatures, the friction coefficient caused 
by a superconducting sample is an order of magnitude 
larger than that of an insulating sample, in contradic- 
tion with the previous theoretical prediction that NCF 
generically scales with the resistivity of the sampleP Fur- 
thermore, a universal feature has been identified in these 
experiments,^ namely, at low temperatures, while the 
induced spring constant increases monotonically with de- 
creasing distance, the friction coefficient displays a max- 
imum at certain distance of a few nanometers. Such fea- 
ture has been consistently observed in insulating materi- 
als as well as superconducting materials, both below and 
above the superconducting transition temperature. 

In this paper, we propose an explanation for these ex- 
perimental findings.^ Specifically, we propose that the 
observed behavior of NCF is due to the relaxation dy- 
namics of the surface defects. 

i) First we lay out the general formalism. It is illus- 
trated by considering localized spins. We assume that 
spins remain unscreened due to the insulating or super- 
conducting gap. We assume that cantilever also carries 
localized spins that couple to the collection of spin sites 
on the surface of the sample which have characteristic re- 
laxation time Td as function of tip-sample distance d. We 
model this relaxation time as a collection of Debye relax- 
ors with distributed relaxation times that results in a typ- 
ical glass-like backaction dynamics, see below. Analysis 
of the experimental data indeed indicates a distribution 
of relaxation timesP^ Alternatively defects can be local 
charged sites. In that case the charging and uncharging 
dynamics of the defects leads to the random electric field 
probed by a cantilever. 

ii) Using our general formalism we then extract from 
experimental data the distance dependence of tip-sample 
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coupling A[d) and relaxation time Td- A{d) can be fitted 
by simple power laws, but Td displays surprising behav- 
ior. For extremely hard cantilever at low temperature, 
clear divergence in Td is observed. Weaker singularities 
are seen for moderately hard cantilevers. We propose an 
explanation for such behavior. 

iii) We show next that the combined effect of the vari- 
ation of the tip-sample coupling and the relaxational dy- 
namics on the sample surface explains the observed dis- 
tance dependence of the friction coefficient and the in- 
duced spring constant. We also estimate the order of 
magnitude of the friction coefficient, and the experimen- 
tally observed values can be easily achieved in our frame- 
work. 

Formalism: We model the tip of the cantilever by a 
massive particle with effective mass m, moving in a one- 
dimensional harmonic potential V = ^kx 2 . Here k is 
the spring constant and x the displacement. We con- 
sider that there are some randomly distributed active de- 
grees of freedom on the sample surface that interact with 
the tip. Generically such interaction will then produce 
a frictional force, damping the motion of the tip. The 
tip motion is thus governed by the generalized Langevin 
equation 
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t')-x(t')dt> = F x (t), (1) 



with 7(f) the backaction term that we write as a dynami- 
cal damping. There is also a residual random force F x (t) , 
which fluctuates rapidly. For an equilibrium system, the 
frictional force and the random force are conn ected by 
the fluctuation-dissipation theorem (FDT)JEEJl 
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where 7(w) is the Fourier transform of j(t), i.e. ,7(0;) = 
J Q °° dte~ l0Jt j(t). j(oj) is generally a complex function. 
One can write it as j(oj) = T(ui)— ^fc; nt (w). The real part 
r describes the effect of dissipation, and the imaginary 
part leads to a modification of the spring constant. T = 
To + Tint, with Tq the intrinsic cantilever friction and T^t 
the NCF resulting from tip-sample interaction. 

Spin-spin interaction: The mechanism we consider is 
quite general. It applies to interactions in different chan- 
nels, e.g. spin or charge. To be concrete, we will study 
in detail the spin-spin interaction. Since essentially the 
same relaxational behavior was observed in superconduc- 
tors (NbSe 2 ) and insulators (SrTiOs), 10 spin-spin inter- 
action is also a plausible choice. We notice that spin 
relaxation of surface magnetic defects is regarded as the 
origin of f/f flux noise in superconducting devicesP^HIH 

We consider there are randomly distributed localized 
spins on the sample surface and they interact with the 
spin localized on the tip, with a Hamiltonian H = 
Y^i a Ji( x )StipSi- Here 5" ; is the spin operator on the 
tip, 5? the spin operator on the sample surface, and the 



coupling J? can be of different types, e.g. Ising, Heisen- 
berg. The force in the x direction is F x = —dH/dx. j(oj) 
thus reads 
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Y,(Jn^J^^ ip (0)SnO)S b tip (t)S^t)). (3) 
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The spins on the sample surface interact with each 
other, and the tip spin provides an external magnetic 
field H for the surface spin system. The Hamiltonian 
of the surface spin system is thus of the form, H sur f = 

- Ei^j J ij s t s j ~ Ei H i S i- Effect of the external field 
on dynamics of spins is that it will change the relax- 
ation time r, see below. The surface spin system displays 
disordered behavior, where the cross correlations vanish 
and the dynamics is characterized by the autocorrelation 
function, with (5?(t)5?(0)) = q(t)5 %r 

We treat the dynamics of the spins on the sample sur- 
face is much faster than that of the tip spin due to the 
significantly higher density of scattering centers at the 
surface, hence cantilever spin can be viewed as static. 
7(0;) can then be factorized into two parts, 



7(0;) 



k B T 



A{d)S d {uj) 



(4) 



The prefactor C s comes from spin degeneracy. The 
frequency-independent tip-sample coupling A(d) = 

lim x _»o(Ej iwx^i( x ^) )' increasing monotonically with 
decreasing d. Introducing a defect density p(r) = 
Ei ^( r — R-i ) , which has average value (p) = ni mp , one 
obtains A(d) = (J drp(r) J'(r) 2 ) = n imp J dr J'(r) 2 . Here 
we have defined J'(r) = £\ lim x ^ -§^Ji{x)5{v - Rj). 

The surface spin susceptibility S d {uj) = 
f °° dte- iut {St (0)5?(i)) is related to the response 

function C{iui n ) = /<f dre'^iSf (0)5? (t)) by the 
classical form of the fluctuation-dissipation theorem 
5R5(w) = (k B T/u)$SC(u), or S(u) = -i(k B T '/w)C(w). 

The interactions among the surface spins are random, 
leading to glassy behavior. Experiments also indicate a 
distribution of relaxation timesP^l We thus assume C(u) 
to have the usual phenomenological form typical for a 
glass system C(ui) = Cq/(1 — iujTd) a , with the exponent 
< a < 1, and the relaxation time T^PiHII] 7(0;) now 
takes the form 

l[uj) = _ i 9l A{d) 1 (5) 

LU (1 - lUJT d ) a 

For cantilevers with high quality factor,^ the intrin- 
sic cantilever friction can be ignored. Defining Aid) = 
C t A(d), with Ci = C s C , the NCF thus reads 

r ( \ M nSin[aarctan(wr d )] 
uTM = A(d) (1 + (w7a)a)o/a , (6) 

and the induced spring constant is 



fci nt (w) = A(d) 



cos [a arctan(wTd)] 

(1 + [uT d ) 2 y/ 2 



(7) 
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FIG. 2: (Color online) Relaxation time, tip-sample coupling and friction at low temperature T = 4.2K for i) NbSe2 and 
fo = 300kHz, ii) NbSe2 and fo — 31.6kHz, iii) SrTiOa and fo = 31.6kHz, from top to bottom. The star and diamond are 
experimental values (data by Saitoh et alP^). The dashed lines denote where Tint has a maximum. We use Eqj6]tofit r int . The 
glass exponent is chosen to be a — 0.9. Since Hd oc y/Ad, we use Ad as a measure of the magnetic field Hd, and define Hd = y/Ad- 

For case i), ojT d = Cl.r \H d - Hsp\ , with Cl = 0.12, Cr = 0.04. For cases ii), iii), ujr d = Cexp( — (H d — Hsp) /Hq) with 

C {u) = 0.6, H]$ = 9,H^ i} = 11.4 and C* ( " l) = 0.1, = 4,H { " l) = 8.9. The tip-sample coupling A(d) = A (d 2 + rg) -a , 

with a« = 3,4° = 2.3 x lOV^ = 4.8; a ( "> = 1.5, A^ = 6.0 x lO 3 ,^" 5 = 2.2; and = 3, A { iH) = 7.5 x 10 7 ,r^ ll) = 8.2. 

We notice that for case i) there is a jump in Ad, which is "inherited" from the singularity in uiTd, and would disappear in a 
more realistic modeling. The high temperature result is included in Appendix II. 



Their ratio 

wr; nt /feint = tan[a arctan(wr [ j)] (8) 

depends only on the relaxation dynamics of surface de- 
fects. We first extract turd from the experimental data 
using EqjHJ and then calculate the tip-sample coupling 
A(d) from Eq(6]or Eq. The results are shown in Fig. 

12 

We approximate the effect of the tip spin as produc- 
ing an uniform magnetic field Hd on the sample surface. 
When the spin-spin interaction decays with distance as 
J (I) ~ l~ a , one has A(d) ~ (d 2 + rg)" a , and H d ~ 
{d 2 +rl) a ' 2 , with r a cutoff. If J (I) ~ exp(-l 2 /r\), then 
A(d) ~ exp{-2d 2 / r 2 A ) and H d ~ exp(-<i 2 /r^). A(d) can 
be fit with simple power laws (see Fig. |2jc)(g)(k)). 

Noticing Hd ~ \J A(d), we also extract the field de- 
pendence of relaxation time in Fig. [2]jb)(f)(j). The most 
surprising result is that for NbSe2 probed by extremely 
hard cantilever (resonance frequency fo = 300kHz), the 
relaxation time shows clear divergence when approach- 
ing certain field strength. This indicates that the sur- 
face spin system falls into the mean field university class 
and is consistent with behavior represented by the long- 
range, weakly interacting Husimi-Temperley model (sec 
Appendix I) . In this model, as one increases the magnetic 
field, the free energy changes from a double well structure 
to a single well structure, and the metastable minimum 



and the barrier combine at certain value of field strength 
to form a saddle point, known as the spinodal point. Near 
the spinodal point H — Hgp , the relaxation of the system 
slows down dramatically, t ~ \H — i?sp| ^ 2 P^^ We 
use these predictions to model spin behavior here. The 
experimental result of lot for NbSe2 with /q = 300kHz 
can be fitted by such mean field form (see Fig. [2^b)). 
In models with shorter range interactions, the change in 
relaxation time is smeared] 21 " 23 ^ as can be seen in Fig. 
[2jF)(j) for the two cases with f = 31.6kHz. They can be 
fitted by an exponential, r = tq exp(—(H — Hgp) 2 / Hq) , 
for H > H SP . 

The maximum in Fj nt can be understood as a result 
of enhanced relaxation time near the spinodal point. In 
the region lot -C 1, while fcj nt ~ A(d) increases monoton- 
ically with decreasing d, I^nt ~ aA(d)ujTd is determined 
by the competition between A{d) and r d - When is sin- 
gular, Tint is determined predominantly by Td near H$p, 
and the maximum is located right at the spinodal point. 
When the singularity in is smeared, due to the dis- 
tance dependence of tip-sample coupling, the maximum 
moves away from the spinodal point towards smaller d 
(see Fig. gd) (h)(1) for numerical fit). 

Let us now estimate the order of magnitude for the 
friction coefficient. Friction is essentially determined by 
the following three factors: the tip-sample interaction 
energy J, the surface impurity density nj mp , and the 
characteristic energy scale of the surface susceptibility 
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E c ~ l/{Sd(ui))- Including spin degeneracy, one can 
write the friction term as uj Ti Dt ~ s 2 S 2 ni mp J 2 /Ec- The 
prefactor s 2 S 2 is about 1-10. Assuming there are about 5 
impurities per square nanometer, to get the experimental 
value of wor^t, which is about lOOiV/m for d ~ 2nm^ 
one needs to have J 2 /Ec — 10 — lOOeU. If we take 
E c to be of order kgT, where T — 4.2 is the tem- 
perature at which the experiment is performed,^ then 
E c ~ 3 x 10~ 4 eU, and the coupling is J ~ 3 — 30meU, 
which can be achieved. 

Charge-charge interaction: For completeness, we also 
consider the possibility that friction arises from interac- 
tions in the charge channel, though this mechanism may 
not apply foi^ . We notice that, as was shown irP, 
coupling between charge on the cantilever tip and ion 
vibrations on the sample surface can produce strong en- 
hancement of NCF. This me chanism may be responsible 
for the observed NCF inPaH 

Essentially same logic outlined earlier applies with ob- 
vious substitution of Coulomb interactions for spin in- 
teractions. We assume there are some randomly dis- 
tributed two-level fluctuators (TLF) on the sample sur- 
face, formed from localized electronic trap statesP^H 
The charge on the tip of the cantilever provides an exter- 
nal electric field, favoring one of the two states. The 
TLF are now governed by the Hamiltonian, T~L S urf — 
~ Ei^j JijQiQj - J2i ViQi, witn charge Qi = 0, 1, the 
coupling Jij a random number, and Vi the electric po- 
tential created by the tip charge. 

The result is qualitatively the same as the case with 
spin-spin interactions, though the order of magnitude can 
be different. Here u>oTi nt ~ q 2 n-i mp J 2 e ^/E^\ For d ~ 
2nm, the Coulomb potential is Ju\ ~ 0.3eU. Taking 

Eq ~ ksT ~ 3 x 10~ 4 eU, to get the experimental result 
of cJoTint — IOOTV/to, one needs to have q 2 n imp ~ 2, 
with q the tip charge in unit of the elementary charge 
e, and n lmp the number of surface charges per square 
nanometer. This can be easily achieved. 

Conclusion: In conclusion, we have proposed a gen- 
eral mechanism to explain the distance dependence of 
the friction coefficient and the induced spring constant 
of an oscillating cantilever. A universal ingredient of the 
proposed mechanism is the backaction effects of relax- 
ational dynamics of the defects on the sample surface. 
This mechanism also explains nicely the observed order 
of magnitude of the friction coefficient. Furthermore, our 
formalism provides a general framework for experimen- 
talists to extract separate information about tip-sample 
coupling and surface dynamics, thus enabling more de- 
tailed investigation of surface properties. One way to test 
our theory would be to examine the magnetic field depen- 
dence of the friction coefficient and the induced spring 
constant. 
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Appendix I: Divergence of relaxation time near the 
spinodal point 

We consider here the long-range, weakly interacting 
Husimi- Temper ley model where each spin <7j = ±1 inter- 
acts equally with every other spin. In the presence of an 
external magnetic field H , the Hamiltonian reads 

J N 

K = -^M 2 -HM, M = J>, (9) 

i=l 

When H = 0, the free energy / has two stable min- 
ima, and this model displays a second-order phase tran- 
sition. In the presence of a weak magnetic field, one 
minimum becomes a metastable state. At some criti- 
cal field strength, this metastable state becomes unsta- 
ble. At this value of field strength, the potential changes 
from a double well structure to a single well structure, 
and the metastable minimum and the barrier combine 
to form a saddle point, known as the spinodal point 
(see Fig{3|. This point is determined by the condition: 
df /dm = 0, d 2 f /dm 2 — 0, with m = M/N the magneti- 
zation per spin. 

Near the spinodal point, the relaxation of the system 
slows down dram aticall y. For H > H^p, the relaxation 
time is of the for m 20 * 28 ! 

r ~ * exp (b \H H SP f 2 ) , (10) 

when A = (3N 2 / 3 \H - H SP \ > 1. Here b = 
(4/3)/3A7 i/j(/3J- 1). Consider the double scaling limit 
N -> oo, \H - H SP \ -> with A large but finite, the 
relaxation time reduces to t ~ \H — Hsp\ ■ For 
H > Hsp . there is a critical divergence of the relaxation 
timeP 

t ~ \H - H SP \~ 1/2 . (11) 

In models with shorter range interactions, the change 
in relaxation time is smeared. 

Appendix II: Friction at room temperature 

Using the same method as in the main text, we extract 
here the relaxation time and tip-sample coupling at high 
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FIG. 3: Evolution of potential energy with magnetic field in 
Husimi-Temperley model. We are plotting f(m) — —^m 2 — 
Hm + i (i±S i n 1±S + l^m m i^n) pa with /3 = f, J = 
1.32, # = 0,0.03,0.11,0.2. 
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FIG. 4: Relaxation time and tip-sample coupling for NbSe2 at 
room temperature and fo = 31.6kHz. The star and diamond 
are experimental values (data by Saitoh et alP^ ). The solid 
line is a fit to A(d) = ao + aid with ao = 43.2, ai = —9.1. 



temperature. The results are shown in Fig|4j We can 
see that both of them display qualitatively different be- 
havior as compared to the low temperature case. In the 
region where the experimental data can be trusted, i.e. 
d < 4nm, ujTd is essentially distance independent; Ad sets 
in at certain value of tip-sample distance ~ 4.5nm, and 
then grows linearly in distance. These two features can 
be understood by simply assuming that at high tempera- 
tures, there exists certain viscous cloud above the sample 
surface, extending up to ~ 4.5nm high. The NCF that's 
at work at low temperatures is now suppressed by tem- 
perature effects, and the observed friction is due to the 
cantilever interacting with such cloud. Since it is always 
the same cloud, the relaxation time should not change 
with distance. Friction is proportional to the length of 
the part of the cantilever that is immersed in this cloud. 



1 B. Bhushan, Nanotnbology and Nanomechanics: An Intro- 
duction (Springer- Verlag, Berlin, 2008). 

2 I. Dorofeyev, H. Fuchs, G. Wenning, and B. Gotsmann, 
Phys. Rev. Lett. 83, 2402 (1999). 

3 B. C. Stipe, H. J. Mamin, T. D. Stowe, T. W. Kenny, and 
D. Rugar, Phys. Rev. Lett. 87, 096801 (2001). 

4 S. Kuehn, R. F. Loring, and J. A. Marohn, Phys. Rev. 
Lett. 96, 156103 (2006). 

5 J. A. Sidles, J. L. Garbini, K. J. Bruland, D. Rugar, 
O. Ziiger, S. Hoen, and C. S. Yannoni, Rev. Mod. Phys. 
67, 249 (1995). 

6 F. J. Giessibl, Rev. Mod. Phys. 75, 949 (2003). 

7 A. I. Volokitin and B. N. J. Persson, Phys. Rev. B 65, 
115419 (2002). 

8 A. A. Chumak, P. W. Milonni, and G. P. Berman, Phys. 
Rev. B 70, 085407 (2004). 

9 A. I. Volokitin and B. N. J. Persson, Rev. Mod. Phys. 79, 
1291 (2007). 

K. Saitoh, K. Hayashi, Y. Shibayama, and K. Shirahama, 
Phys. Rev. Lett. 105, 236103 (2010). 

1 R. Kubo, M. Toda, and N. Hashitsume, Statistical Physics 
II: Nonequilibrium Statistical Mechanics (Springer- Verlag, 
1985). 

2 R. Kubo, Rep. Prog. Phys. 29, 255 (1966). 

3 R. H. Koch, D. P. DiVincenzo, and J. Clarke, Phys. Rev. 
Lett. 98, 267003 (2007). 

4 R. de Sousa, Phys. Rev. B 76, 245306 (2007). 



15 S. Sendelbach, D. Hover, A. Kittel, M. Muck, J. M. Mar- 
tinis, and R. McDermott, Phys. Rev. Lett. 100, 227006 
(2008). 

16 K. S. Cole and R. H. Cole, J. Chem. Phys. 9, 341 (1941). 

17 D. W. Davidson and R. H. Cole, J. Chem. Phys. 18, 1417 

(1950) . 

18 D. W. Davidson and R. H. Cole, J. Chem. Phys. 19, 1484 

(1951) . 

19 K. Binder, Phys. Rev. B 8, 3423 (1973). 

20 T. Mori, S. Miyashita, and P. A. Rikvold, Phys. Rev. E 
81, 011135 (2010). 

21 J. E. McDonald, Am. J. Phys. 30, 870 (1962). 

22 J. E. McDonald, Am. J. Phys. 31, 31 (1963). 

23 P. A. Rikvold, H. Tomita, S. Miyashita, and S. W. Sides, 
Phys. Rev. E 49, 5080 (1994). 

24 A. I. Volokitin and B. N. J. Persson, Phys. Rev. Lett. 94, 
086104 (2005). 

25 M. Kisiel, E. Gnecco, U. Gysin, L. Marot, S. Rast, and 
E. Meyer, Nature Materials 10, 119 (2011). 

26 C. T. Rogers and R. A. Buhrman, Phys. Rev. Lett. 55, 
859 (1985). 

27 M. Constantin and C. C. Yu, Phys. Rev. Lett. 99, 207001 
(2007). 

28 W. Paul, D. W. Heermann, and K. Binder, J. Phys. A 22, 
3325 (1989). 



